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Abstract 

A  two-dimensional  numerical  model  has  been  established  to  investigate  the  performance  of  the  PEM  fuel  cells.  Parameters  used  in  the 
analysis  include  the  porosity  and  thickness  of  the  gas  diffuser  layer  (GDL).  Results  show  that  increasing  the  porosity  of  gas  diffusion  layer 
causes  the  increasing  of  mass  transfer  of  fuel  and  air  and  results  in  a  higher  reaction  rate.  Therefore,  a  better  performance  of  the  fuel  cell 
and  more  fuel  consumption  rate  are  observed.  It  is  also  demonstrated  that  the  performance  of  the  fuel  cell  increases  with  a  decrease  in  the 
thickness  of  gas  diffusion  layer.  The  effects  of  liquid  water  condensation  and  flow  directions  of  fuel  and  air  are  also  considered  in  this  analysis. 
Predicted  results  show  that  the  performance  of  the  PEM  fuel  cell  without  consideration  of  liquid  water  effect  is  always  higher  than  that  with 
consideration  of  liquid  water  effect.  In  addition,  the  performance  of  fuel  cell  with  co-flow  pattern  of  fuel  and  air  is  larger  than  that  with  counter 
flow. 

©  2005  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  is 
considered  as  a  most  promising  and  environmental  friendly 
alternative  power  source  for  both  stationary  and  mobile  appli¬ 
cations  due  to  its  high  efficiency  and  low  operating  temper¬ 
ature.  Conventional  PEMFCs  operate  on  hydrogen  as  a  fuel 
and  oxygen  or  air  as  an  oxidant  to  produce  electricity,  and 
heat  and  water  are  the  major  by-products.  The  gas  streams 
have  to  be  humidified  to  keep  the  membrane  water  swollen 
in  order  to  obtain  sufficiently  high  protonic  conductivity.  The 
reactant  gases  in  the  channel  transport  through  the  gas  dif¬ 
fuser  layer  (GDL)  and  reach  the  catalyst  layer  to  carry  out 
electrochemical  reaction.  In  past  decades,  substantial  efforts 
have  been  focused  on  fuel  cell  structure  design,  development 
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of  catalyst,  structure  of  catalyst  and  gas  diffuser  layer  and  the 
development  of  membrane  with  high  performance,  thermal 
and  water  management  [1-8]. 

The  performance  of  the  fuel  cell  depends  on  the  kinetics  of 
the  electrochemical  process  and  performance  of  the  compo¬ 
nents.  Many  experiments  have  been  conducted  by  researchers 
for  varying  operating  parameters.  However,  experimental 
investigations  are  costly.  There  are  several  models  explain 
the  mass  transfer  effect  for  different  species  in  the  fuel  cell 
developed.  The  models  can  predict  the  cell  performance  pre¬ 
cisely.  Several  mathematical  modeling  approaches  have  been 
proposed  to  describe  the  transport  phenomena  within  PEM¬ 
FCs  and  the  impact  of  the  water  condensation  on  the  cell 
performance.  In  recent  year,  two-dimensional  models  [9-12] 
were  developed  in  order  to  investigate  the  distribution  of  gas 
components,  water  content  of  membrane,  current  density  and 
net  water  transport  along  the  flow  channel.  Gurau  et  al.  [13] 
presented  a  comprehensive  model  for  the  entire  sandwich 
of  a  PEMFC  including  the  gas  channels  and  considered  the 
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Nomenclature 

Ajo  reference  exchange  current  density  (Am-2) 
a  chemical  activity  of  water  vapor  in  cathode 

Cp  quadratic  drag  factor 

C  concentration 

D  diffusivity  (m2  s-1) 

F  Faraday  constant,  96485  c  mol-1 
i  current  density  (A  m-2) 

j  current  density  (A  m-2) 

k  permeability  (m2) 

kc  condensation  rate  constants 

&e  evaporation  rate  constants 

M  molecular  weight 

P  pressure  (atm) 

Pi  partial  pressure  for  i  species  (atm) 

R  universal  gas  constant  (8.314  J  mol-1  K-1) 

R\ s  volume- to- surface  ratio  of  the  porous  material 

s  the  ratio  of  the  volume  of  pore  occupied  by 

liquid  water  to  the  volume  of  pore  in  the  porous 
medium 

S  source  term  in  momentum  equation 

Sc  source  term  of  chemical  reaction  in  the  species 

concentration  equation 

Sj  source  term  in  phase  potential  equation 

Sl  source  term  with  consideration  of  liquid  water 
in  the  species  concentration  equation 
t  thickness  (m) 

T  temperature  (K) 

U,  V  velocities  in  the  X  and  Y direction  (ms-1) 

X,  Y  rectangular  coordinate  system  (m) 

Z  number  of  electrons  transferred 

Zf  charge  transfer  coefficient 

Greek  letters 

a  charge  transfer  rate 

£  porosity 

0  membrane  potential 

v  kinematic  viscosity  (m2  s-1) 

a  electric  conductivity  (Q ~ 1  m- 1 ) 

p  density  (kg  m- 3 ) 

r  tortuosity  of  the  pore  in  the  porous  medium 

0  overpotential 

Superscripts 

ref  reference  value 

Subscripts 

a  quantity  in  anode 

c  quantity  in  cathode 

eff  effective  value 

g  of  gas  diffuser  layer 

H+  for  proton 

H2  for  hydrogen 


H2O  for  water 

i  for  i  species 

m  of  membrane 

O2  for  oxygen 

sat  saturation  pressure  for  water  vapor 

total  for  total  value 

v  in  the  X  direction 

y  in  the  Y  direction 


gas-liquid  phases  in  separate  computation  domains  for  trans¬ 
port  in  the  gas  distribution  channels.  It  was  assumed  that 
the  catalyst  layer  is  infinitesimally  thin.  Hsing  and  Futerko 
[14]  developed  a  two-dimensional  model  of  coupled  fluid 
flow,  mass  transport  and  electrochemistry  of  a  PEMFC  with 
taking  into  account  the  dependence  of  diffusion  coefficient 
of  liquid  water  in  membrane.  However,  their  model  does 
not  resolve  catalyst  layers  and  hence  ignores  the  influence 
of  spatial  nonuniformity  of  water  content  on  catalyst  layer 
performance.  Yi  and  Nguyen  [15]  proposed  an  along-the- 
channel  model  for  evaluating  the  effects  of  various  design  and 
parameters  on  the  performance  of  a  PEMFC.  The  results  show 
that  the  humidification  of  the  anode  gas  is  required  to  enhance 
the  conductivity  of  the  membrane,  and  the  liquid  injection 
and  higher  humidification  temperature  can  improve  the  cell 
performance  by  introducing  more  water  into  the  anode.  It 
is  noted  that  the  mass  transport  processes  in  the  presence 
of  liquid  water  are  not  considered  in  these  two-dimensional 
models. 

In  a  typical  PEMFC,  the  gas  diffuser  layers  (GDLs)  are 
required  to  have  both  reactant  gases  supply  to  the  catalyst 
layer  and  water  supply  and  removal  in  either  vapor  or  liq¬ 
uid  form.  Although  it  is  a  seemingly  minor  component  in 
a  fuel  cell,  GDL  is  one  of  the  important  parts  of  the  PEM 
fuel  cell.  Any  change  in  the  composition  or  the  morphology 
of  the  GDL  can  lead  to  a  substantial  influence  on  fuel  cell 
performance.  Jordan  et  al.  [16]  experimentally  examined  the 
influence  of  diffusion-layer  morphology  on  cell  performance. 
They  developed  a  model  of  the  hydrophobilicity  and  porosity 
of  the  diffusion  layer  to  explain  the  influence  of  the  diffusion- 
layer  morphology.  Later,  they  [17]  obtained  results  which 
showed  that  a  low-porosity  acetylene  black  enhances  water 
removal  from  the  membrane  electrode  assembly  (ME A),  and 
thus  gives  an  improvement  in  gas-diffusion.  Recently,  Lee 
et  al.  [18]  also  studied  experimentally  the  effect  of  the  fab¬ 
rication  method  and  the  thickness  of  the  GDL  on  the  cell 
performance.  However,  their  focus  was  on  the  fabrication 
method,  rolling,  spraying  and  screen  printing,  of  GDL  effect¬ 
ing  on  the  cell  performance.  Gurau  et  al.  [19]  proposed  a 
mathematical  model  to  obtain  an  analytic  solution  of  the  mass 
transport  of  reactant  gas  in  a  half-cell,  in  which  the  effects 
of  the  porosity  and  the  tortuosity  of  the  GDL  and  catalyst 
layer  were  explored  due  to  the  fact  that  the  pores  may  be 
partially  filed  with  liquid  water.  In  their  findings,  the  cata¬ 
lyst  layer  resistance  at  low  current  densities  is  equal  to  half 
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the  value  of  the  catalyst  layer  resistance  corresponding  to  the 
limiting  current  density.  Chu  et  al.  [20]  also  used  a  half-cell 
model  to  investigate  the  effect  of  the  change  of  the  GDL  on 
the  performance  of  a  PEMFC.  They  found  that  a  fuel  cell 
embedded  with  a  GDL  with  a  larger  averaged  porosity  will 
consume  a  greater  amount  of  oxygen,  so  that  a  higher  cur¬ 
rent  density  is  generated  and  a  better  fuel  cell  performance  is 
obtained.  Natarajan  and  Nguyen  [21]  developed  a  transient, 
two-dimensional  mathematical  model  to  study  the  effect  of 
thickness  of  GDL  for  the  cathode  of  a  PEMLC.  However,  their 
model  does  not  resolve  the  catalyst  layer  and  water  transport 
in  the  membrane  is  not  considered.  They  found  that  the  per¬ 
formance  of  the  cathode  was  dominated  by  the  dynamics  of 
liquid  water,  especially  in  the  high  current  density  range,  and 
a  better  performance  can  be  observed  when  a  thinner  GDL  is 
utilized  in  the  fuel  cell. 

In  the  study  of  the  PEMFCs,  the  water  management  is 
one  of  the  most  important  issues.  Nguyen  and  White  [22] 
proposed  a  model  of  the  water  and  heat  management  of 
the  PEMFC  systems,  which  includes  the  effect  of  electro¬ 
osmosis,  diffusion  of  water,  heat  transfer  from  solid  phase 
to  gas  phase  and  latent  heat  as  water  evaporation  and  con¬ 
densation.  They  found  that  the  ohmic  loss  is  considerable 
at  high  current  and  the  voltage  loss  is  twice  amount  of 
that  of  the  cathode  electrode.  Baschuk  and  Li  [23]  devel¬ 
oped  a  mathematical  model  with  variable  degrees  of  water 
flooding  in  the  PEMFC.  Physical  and  electrochemical  pro¬ 
cesses  occurring  in  the  membrane  electrolyte,  the  cathode 
catalyst  layer,  the  electrode  backing  layer  and  the  flow  chan¬ 
nel  were  considered.  They  discovered  that  when  the  air  was 
used  as  the  cathode  fuel,  the  flooding  phenomena  are  sim¬ 
ilar  at  different  operating  conditions  of  the  pressure  and 
temperature  comparing  with  experimental  results.  Wang  et 
al.  [24]  have  numerically  modeled  the  two-phase  flow  and 
transport  in  the  air  cathode  of  PEMFCs.  They  found  that 
the  transport  processes  of  liquid  water  and  water  vapor  are 
affected  by  capillary  effect  and  molecular  diffusion,  respec¬ 
tively.  A  two-dimension  model  was  developed  by  Ge  and 
Yi  [25]  to  investigate  the  effects  of  operation  conditions 
and  membrane  thickness  on  the  water  transport.  In  their 
study,  the  liquid  water  effect  on  the  effective  porosity  for 
gas  transport  was  considered  to  simplify  the  model  of  the 
two-phase  flows  in  porous  layers.  The  results  revealed  that 
the  cell  performance  can  be  enhanced  by  increasing  the  cell 
temperature. 

From  the  literatures  reviews  cited  above,  it  is  found  that  a 
full-cell  mathematical  model  with  consideration  of  catalyst 
layer  and  membrane  water  transport  in  a  PEMFC  has  not 
been  well  examined  yet.  This  motivates  the  present  study. 
The  objective  of  this  work  is  to  establish  a  two-dimensional, 
full-cell  mathematical  model  with  consideration  of  water  sat¬ 
uration  to  investigate  the  effects  of  both  the  porosity  and 
dimension  of  gas  diffuser  layer  on  the  performance  and  gas 
transport  phenomena  of  a  PEMFC.  Additionally,  the  effects 
of  flow  direction  on  the  cell  performance  are  taken  into 
account  in  the  analysis. 


anode  MEA  cathode 


Fig.  1.  Schematic  diagram  of  complete  PEMFC  and  coordinate  system. 

2.  Analysis 

The  PEM  fuel  cell  model  described  in  this  study  is  devel¬ 
oped  to  investigate  the  effects  of  porosity  and  thickness  of 
GDL.  The  computational  domain  includes  a  membrane  sand¬ 
wiched  between  two  gas  diffusion  electrodes,  and  the  flow 
channels  of  both  anode  and  cathode.  A  schematic  illustration 
of  a  PEMFC  in  the  co-flow  mode  and  associated  coordinate 
system  is  shown  in  Fig.  1.  U  and  V  are  the  velocity  compo¬ 
nents  in  the  X  and  Y  directions,  respectively.  The  PEMFC  in 
the  counter-flow  mode  is  also  analyzed  to  study  the  perfor¬ 
mance  of  a  PEMFC  between  co-flow  and  counter-flow  modes 
for  a  comparison.  Consequently,  to  simplify  the  problem,  a 
steady  state,  two-dimensional,  multi-species  and  along-the- 
channel  model  of  a  full-cell  PEMFC  is  employed  for  the 
analysis.  Four  species  were  considered:  hydrogen,  oxygen, 
nitrogen  and  water  vapor.  Stationary  conditions  is  assumed 
in  this  fuel  cell,  also  the  effect  of  gravity  is  neglected.  It  is 
assumed  that  the  electrochemical  reactions  takes  place  only 
in  the  catalyst  layer,  and  the  gas  mixtures  in  the  flow  chan¬ 
nels  are  also  considered  to  be  perfect  gases.  Based  on  the 
definition  of  the  Reynolds  number  and  the  velocity  used  in 
this  work,  the  flow  in  the  fuel  cell  is  laminar.  Therefore,  all 
the  transport  equations  were  formulated  for  laminar  behav¬ 
ior.  The  GDL,  catalyst  layer  and  PEM  are  assumed  to  be 
isotropically  porous  materials. 

In  the  development  of  this  numerical  model  of  the  PEM 
fuel  cell  as  well  as  corresponding  analyses,  are  crucial  to 
gain  a  good  understanding  of  the  effect  of  GDL  porosity  and 
thickness  on  the  cell  potential.  Most  of  the  previous  studies  in 
this  area  have  assumed  for  simplicity  constant  porosity  of  the 
GDL.  This,  however,  does  not  reflect  to  the  physical  reality. 
Although  it  has  been  recognized  that  the  performance  of  the 
fuel  cell  can  be  significantly  influenced  due  to  the  porosity 
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variation,  systematical  studies  are  still  rarely  reported  in  the 
literature.  In  this  study,  porosity  of  GDL  is  chosen  as  0.3, 
0.4,  0.5  and  0.6.  These  values  are  chosen  as  closed  as  to 
the  realistic  situation  and  to  reach  a  quantitative  influence  of 
GDL  porosity  on  the  cell  performance. 


Species  equation: 


£eff 


-\-Sc  + 
(9) 


2.7.  Governing  transport  equations 


According  to  the  descriptions  and  assumption  above,  the 
basic  transport  equations  for  the  two-dimensional  PEM  fuel 
cell  are  given  as  the  following: 


Continuity  equation: 


dU  dV 
—  +  —  =0 
dx  dY 

Momentum  equation: 

du 
u 


u 


dx 

dv 

dx 


du 

1  dp  , 

fd2U  d2U 

+  v 

dY 

T  v 

pdX  ' 

\ dx 2  +  dY2 

dV 

1  dP  , 

2d2V  d2V 

T  V  — 

dY 

T  v  ( 
p  dY  \ 

v  dX 2  +  3  Y2 

+  s 


X 


+  s 


y 


(1) 


(2) 


(3) 


The  source  terms  in  Eqs.  (2)  and  (3)  are  different  in  each 
computational  domain  due  to  the  pressure  difference  when 
fluids  pass  through  porous  medium.  So,  the  Sx  and  Sy  in  the 
GDL  and  the  catalyst  layer  are: 


where  D/?eff  is  the  effective  diffusion  coefficient  of  the  i- th 
composition  of  fuel  reactant.  The  effective  diffusion  coef¬ 
ficient  is  related  to  the  diffusion  coefficient  in  a  nonporous 
system  D;  by  [27]: 

A,e  ff  =  A^eff  (10) 


where  r  is  the  tortuosity  of  the  pore  in  the  porous  medium, 
which  is  1.5  for  GDL  and  catalyst  layer,  and  6  for  the  mem¬ 
brane.  Sc  is  the  source  term  of  chemical  reaction  in  the 
species  concentration  equation,  and  is  0  in  the  gas  flow  chan¬ 
nel  and  GDL.  Lurthermore,  Sc  is  different  with  the  reactant 
gases  in  the  catalyst  layer  and  PEM,  e.g.,  Sc  is  —  for 

hydrogen  in  catalyst  layer,  —  4Fq  { 1  for  oxygen  in  catalyst 


layer,  and 


J  C 


2FC[  otal 


,c 


for  water  vapor  in  catalyst  layer,  and 


w D<,ef/H+CH+  (ff  +  ff  )  in  the  PEM.  The  parameters 

7a  and  jc  indicate  the  current  density  at  the  anode  and  cath¬ 
ode  sides,  respectively,  and  can  be  described  by  the  following 
Butler- Volmer  equations: 


Sx  =  - 


Sy  = 


VSeff 

k 

V£eff 


u  — 


V  - 


c‘’Cr/u  s/WW* 

Xk 

(4) 

h  =  Aioef  | 

( Cn2 

i  >-ref 

Vch2 

e‘"Cr/v 

Xk 

(5) 

jc  =  Aioef  1 

(  Co2 

i  rref 

Vco2 

>Oa  F/RT)  _ 


1 


facF/RT)r 


faa  F/RT)  _ 


1 


JacF/RT)r 


(11) 

(12) 


and  the  Sx  and  Sv  in  the  PEM  are: 


•y 


4  = 

k  V£ 

k  du 

-\ - Zf  F  •  V  0 


V£eff 


dx 


k  Vk 


-\- - Zf  Cy\+  F  •  X  0 

V£eff 


dV 

~dY 


where  £eff  is  the  effective  porosity,  Cp  the  quadratic  drag  fac¬ 
tor  and  k  is  the  permeability,  Zf  the  charge  transfer  coefficient, 
CH+  the  concentration  of  proton,  F  the  Laraday  constant  and 
0  is  the  membrane  potential.  In  the  analysis,  Blake-Kozeny 
equation  [26]  is  used  to  model  k  as  below: 


where  AjQ6f  is  the  exchange  current  density,  aa  and  ac  the 
electric  charge  transport  rates  in  anode  and  cathode  catalyst 
layers,  rj  the  overpotential,  R  the  gas  constant  and  T  is  the 
temperature  of  the  fuel  cell. 

To  consider  liquid  water  effect  in  this  work,  a  simplified 
two-phase  model  is  used  to  describe  the  water  transport  in  the 
PEM  fuel  cell.  In  this  work,  the  liquid  water  effect  is  taken 
into  account  by  modifying  the  mass  diffusivity  due  to  the  liq¬ 
uid  water  filling  the  pores  in  the  porous  media  and  the  liquid 
water  generation  in  the  species  concentration  equation.  When 
the  partial  pressure  of  water  vapor  is  greater  than  the  satura¬ 
tion  pressure  of  water  vapor,  the  water  vapor  is  assumed  to 
condense  and  fill  the  pore  in  the  porous  media.  The  source 
term,  Sl  in  the  species  concentration  equation,  Eq.  (9),  rep¬ 
resenting  the  quality  of  liquid  water  can  be  evaluated  [28]: 


^H20 kc  gff  Jl2°  ( Fh20  ~  Aat),  if  At20  >  Aat 

pRT 

keseffs(PSSit  -  Ph2o),  if  At2o  <  Aat 

(13) 


and  Dip  =  67?vs,  and  R\s  is  the  volume-to-surface  ratio  of  the 
porous  material. 


where  the  M  is  the  molecular  weight  and  the  kc  and  kQ  are  the 
condensation  and  evaporation  rate  constants.  The  saturation 
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pressure  of  water  can  be  expressed  as: 

r>  in-2.1794+0.02953r-9.1837xl0_5r2+1.4454xl0_7r3 
~sat  — 

(14) 

It  should  be  noted  that  the  assumption  of  isothermal  condi¬ 
tion  is  utilized  in  this  analysis  to  simplify  the  model.  However, 
the  saturation  of  water  is  greatly  affected  by  the  temperature 
in  the  fuel  cell,  such  as  the  saturation  pressure  in  Eq.  (14). 
Nevertheless,  when  a  single  cell  is  tested,  the  cell  temperature 
is  usually  controlled  using  heating  rods  and  thermocouples. 
As  a  result,  the  operation  temperature  of  a  single  fuel  cell 
can  be  kept  almost  constant  in  all  the  regions  of  the  working 
fuel  cell.  The  temperature  difference  is  small  enough  during 
cell  operation  that  the  assumption  of  isothermal  condition  is 
appropriate  for  the  analysis. 

In  addition,  the  saturation,  s,  is  defined  as  the  ratio  of  the 
volume  of  pore  occupied  by  liquid  water  to  the  volume  of 
pore  in  the  porous  medium,  then  the  effective  porosity  of 
porous  media  is  modified  to  account  the  liquid  water  effect: 

£eff  =  e(l-*)  (15) 


In  order  to  evaluate  the  distributions  of  the  local  current  den¬ 
sity,  the  phase  potential  equation  should  be  solved: 


d  f  d  f  _ 

dx  vm9x )  +  dY  vm 5y )~S] 


(16) 


where  Sj  =  — ja  in  the  anode,  —jc  in  the  cathode  and  0  in 
the  membrane;  crm  is  the  electric  conductivity  of  the  mem¬ 
brane  which  can  be  calculated  by  the  equation  developed  by 
Springer  et  al.  [29]: 


am{T)  =  cr^f  exp 


and  the  reference  electric  conductivity  is: 
a!f  =  0.0051397.  -  0.00326 


(17) 


(18) 


X  = 


0.043  +  17.81a  —  39.85a2  +  36.0a3,  0  <  a  <  1 


14+  1.4(a  —  1), 


1  <  a  <  3 

(19) 


In  Eq.  (19),  the  a  is  the  ability  of  water  vapor  at  the  cath¬ 
ode  side.  Using  the  following  relations  between  the  phase 
potential  0  and  current  density  /: 


I'X  — 


ly  -  CT] 


d& 

dx 

d& 

3 y 


(20) 


(21) 


Eq.  (16)  can  then  be  reduced  to  be: 


dix  diy  _  . 

3 X+  dY  ~  ■7a’ 


at  anode 


(22) 


di x  di y 

- 1 - -  =  /c ,  at  cathode 

dX  dY  J 


2.2.  Boundary  conditions 


(23) 


Boundary  conditions  for  the  dependent  variables  of  the 
transport  equations  at  the  interfaces  between  different  layers 
of  the  same  domain  are  not  required.  At  the  gas  channel  outlet, 
the  fully  developed  flow  conditions  are  assumed: 


dv  dCi 

U  =  —  =  — -  =0 
dY  dY 


(24) 


The  boundary  conditions  at  the  gas  flow  channel  walls  are: 


dX 


(25) 


In  practical  situation,  the  physical  properties  and  their  gradi¬ 
ents  are  continuous  on  the  interface.  So  the  natural  boundary 
conditions  on  the  interface  are  the  same  velocity,  same  con¬ 
centration  and  the  same  gradients.  At  the  interfaces  between 
the  gas  diffusers  and  the  gas  channels,  the  following  boundary 
conditions  are  used: 


dv 

eeff,x+  yy 


x=x+ 


dv 

dx 


X=X- 


yz=z+  =  yz=z-  (26) 


£&ff,X+ 


dQ 

dX 


x=x+ 


dCi 

dX 


Ci,X=X+  ~  Ci,X=X- 


X=X 


1 


(27) 


The  similar  conditions  are  employed  for  the  interfaces 
between  the  gas  diffusers  and  the  catalyst  layers  and  the 
interfaces  between  the  catalyst  layers  and  membrane  can  be 
expressed  as  follows: 


£eff,X+ 


dV 

dx 


x=x+ 


—  £eff,X- 


dv 

dx 


X=X- 


Vx=x+  —  Vx=x- 

(28) 


£Qff,X+ 


dQ 


dX 

CijX=x+  =  CiX=x~ 


dCi 

—  £eff,X- 

X=X+  °X 


x=x- 


(29) 


The  boundary  conditions  for  the  phase  potential  at  the  inter¬ 
face  between  the  catalyst  layer  and  the  membrane  are  0  =  0  at 
the  anode  side,  and  =  0  at  the  cathode  side.  Because  the 
phase  potential  is  a  linear  distribution  in  the  membrane,  the 
phase  potential  boundary  condition  can  be  written  as  ||  =  0. 


3.  Numerical  method 

The  solution  to  the  governing  equations  is  performed  using 
a  finite  volume  scheme  by  dividing  the  model  domain  into 
a  number  of  cells  as  control  volumes.  In  the  finite  volume 
method,  the  governing  equations  are  numerically  integrated 
over  each  of  these  computational  cells  or  control  volumes. 
The  finite  volume  method  exploits  a  collocated  cell-centered 
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Fig.  2.  Comparison  of  predictions  for  various  grid  systems. 


variable  arrangement,  that  implies  all  the  dependent  variables 
and  material  properties  are  stored  at  the  cell  center.  The  aver¬ 
age  value  of  any  quantity  within  a  control  volume  is  given  by 
its  value  at  the  cell  center. 

The  governing  equations  can  be  expressed  in  the  form  of 
a  generalized  transport  equation: 

v  •  (pUcp  -  r0v0)  =  (30) 

where  0  denotes  the  general  dependent  variable,  r $  the 
exchange  coefficient,  the  source  term,  u  the  velocity  vector 
and  p  is  the  density.  With  the  discritization  of  the  govern¬ 
ing  equations,  the  coupled  finite-difference  equations  can  be 
expressed  in  the  form  of: 

<2p0p  =  OE0E  +  «W0W  +  ^N0N  +  «S0S  +  S(p  (31) 

where  0p  is  the  value  of  0  at  the  current  point  P,  0e.  .  .0s 
stand  for  the  values  of  the  grid  points  adjacent  to  the  point  P 
and  <2p. .  .as  are  known  as  the  link  coefficients. 

In  this  work,  the  non-uniform  grid  system  of  85  x  70  is 
employed  for  the  analysis.  In  order  to  examine  the  grid  inde¬ 
pendence  of  the  predictions,  coarse  and  fine  grid  systems  are 
considered  in  the  preliminary  tests.  Effects  of  the  grid  num¬ 
ber  on  the  predictions  of  local  current  density  are  shown  in 
Fig.  2.  The  maximum  deviations  among  the  computations  on 
the  grids  of  50  x  35,  85  x  70  and  130  x  140  are  less  than  3%. 
Therefore,  the  grid  system  of  85  x  70  points  seems  to  be  suf¬ 
ficient  to  resolve  the  behaviors  of  the  reactant  gas  transport 
in  the  present  PEMFC  model. 

4.  Results  and  discussions 

In  this  work,  the  inlet  conditions  for  anode  and  cathode  are 
that  the  velocity  is  1  m  s-1,  the  temperature  is  333.15  K,  and 
the  relative  humidity  is  100%.  The  total  length  of  the  flow 
channel  is  14  cm,  and  the  cross-section  of  the  flow  channel  is 
1  mm  x  1  mm.  The  thickness  of  the  catalyst  layer  is  taken  to 


Table  1 


Values  of  the  major  parameters 


tg  (m) 

£g 

tm  (m) 

Flow  direction 

1 

0.0004 

0.3 

0.00175 

Co-flow  and  counter-flow 

2 

0.0004 

0.4 

0.00175 

Co-flow  and  counter-flow 

3 

0.0004 

0.5 

0.00175 

Co-flow  and  counter-flow 

4 

0.0004 

0.6 

0.00175 

Co-flow  and  counter-flow 

5 

0.0002 

0.4 

0.00175 

Co-flow 

6 

0.0006 

0.4 

0.00175 

Co-flow 

be  0.0000287  m.  Results  are  obtained  for  cases  presented  in 
Table  1,  where  case  2  is  the  typical  case  for  comparison. 

4. 1.  Effect  of  GDL  porosity 

In  order  to  study  the  effect  of  gas  diffuser  layer  (GDL) 
porosity  on  the  cell  performance,  Fig.  3  shows  the  polariza¬ 
tion  (J-V)  curves  with  various  GDL  porosity  of  a  fuel  cell. 
For  comparison,  the  results  without  liquid  water  effects  in 
the  model  are  also  presented.  An  overall  inspection  on  Fig.  3 
discloses  that  the  effect  of  the  GDL  porosity  on  the  cell  perfor¬ 
mance  is  negligibly  small  at  the  condition  of  high  operating 
voltages.  However,  under  low  operating  voltage  conditions, 
the  GDL  porosity  shows  a  significant  influence  on  the  I-V 
curves.  In  addition,  the  cell  performance  is  increased  with  an 
increase  in  the  GDL  porosity  with/without  consideration  of 
liquid  water  effect.  This  is  because  more  reactant  gases  trans¬ 
fer  into  the  catalyst  layer  with  a  higher  GDL  porosity,  which 
in  turns,  more  chemical  reaction  occurs  and  more  reactant 
gases  consume,  which  results  in  a  better  cell  performance.  It 
is  also  observed  that  the  cell  performance  is  over-predicted 
when  the  liquid  water  effect  is  not  taken  into  account  in  the 
modeling.  This  maybe  due  to  the  fact  that  the  void  in  GDL  is 
filled  with  the  liquid  water,  which  in  turn  causes  the  reduction 
of  mass  transfer.  It  is  also  noted  in  Fig.  3  that  the  deviations 
in  the  cell  performance  between  the  results  with  and  without 
consideration  of  liquid  water  effects  are  small  at  high  operat- 


Fig.  3.  Effects  of  GDL  porosity  on  the  polarization  curves  of  the  PEMFC 
with/without  liquid  water  effects. 
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Fig.  4.  Effects  of  GDL  porosity  on  the  local  current  density. 


in g  voltage  conditions.  This  indicates  that  the  fuel  transport 
in  the  PEMFC  can  be  treated  as  gas  phase  at  high  voltage  con¬ 
ditions.  However,  at  low  voltage  conditions,  the  liquid  water 
effects  on  the  cell  performance  are  remarkable  and  cannot 
be  neglected  in  the  modeling.  This  confirms  the  fact  that  the 
mass  transports  at  lower  operating  conditions  are  significant 
and,  in  turn,  more  water  is  generated  in  the  catalyst  layer  of 
the  cathode  side.  Therefore,  two-phase  flow  effects  should  be 
considered  under  low  operating  voltage  conditions. 

It  is  shown  from  above  that  the  GDL  porosity  has  a 
remarkable  impact  on  the  cell  performance.  It  is  necessary 
to  investigate  the  GDL  porosity  effect  on  the  local  transport 
characteristics  of  the  PEMLCs.  The  GDL  porosity  effects  on 
the  local  current  density  distributions  at  V  =  0.7  and  0.2  V  are 
presented  in  Fig.  4.  To  investigate  the  liquid  water  effects,  the 
results  with  and  without  considerations  of  liquid  water  in  the 
modeling  are  both  presented  in  Fig.  4.  The  local  current  den¬ 
sity  decreases  with  the  axial  direction  and  is  kept  at  almost 
the  same  low  value  under  a  higher  operating  voltage.  It  is 
expected  that  for  the  lower  operating  voltage,  the  electro¬ 
chemical  reaction  is  stronger  in  the  catalyst  layer.  Therefore, 
the  fuel  is  consumed  in  a  higher  rate  when  the  fuel  flows 
along  the  main  flow  direction,  which  in  turn,  causes  the 
substantial  variations  in  the  local  current  density.  A  higher 
local  current  density  is  noted  for  the  PEM  fuel  cell  with  a 
higher  GDL  porosity  near  the  entrance.  This  indicates  that 
an  increase  in  GDL  porosity  enhances  the  diffusion  transport 
of  the  reactant  gas  through  the  porous  layers  and,  in  turn, 
increases  the  local  current  density.  However,  at  the  down¬ 
stream  region,  the  current  density  decreases  as  the  GDL 
porosity  increases.  The  degradation  in  local  current  density 
in  the  downstream  region  is  a  reflection  of  the  more  effi¬ 
cient  fuel  transport  and  the  chemical  reaction  in  the  upstream. 
Besides,  the  predictions  reveal  that  the  consideration  of  the 
liquid  water  in  the  simulation  leads  to  the  relatively  lower 
current  density  than  those  in  the  cases  without  liquid  water 
effect. 


Fig.  5.  Effects  of  GDL  thickness  on  the  polarization  curves  of  the  PEMFC 
with/without  liquid  water  effects. 

4.2.  Effect  of  GDL  thickness 

The  effect  of  the  gas  diffuser  layer  (GDL)  thickness  on 
the  polarization  curves  of  a  PEM  fuel  cell  is  shown  in  Fig.  5. 
It  is  obvious  that  the  cell  performance  with/without  liquid 
water  effect  increases  as  the  GDL  thickness  decreases,  espe¬ 
cially  at  lower  operating  voltage  conditions.  This  is  due  to 
the  increase  of  concentration  gradient  with  a  thinner  GDL, 
which  in  turn,  results  in  a  higher  mass  transfer  to  the  cat¬ 
alyst  layer.  Therefore,  a  larger  current  density  takes  place. 
However,  it  does  not  show  the  same  trend  at  higher  operating 
conditions  ( V >  0.6  V).  Since  the  effects  of  the  GDL  thickness 
show  a  significant  impact  on  the  cell  performance,  the  effect 
of  GDL  thickness  on  the  local  transport  characteristic  is  also 
studied.  Fig.  6  depicts  the  effect  of  the  GDL  thickness  on 
the  local  current  density  along  the  axial  direction  at  V=0.7 
and  0.2  V.  It  is  clearly  seen  that  the  local  current  density 
for  V=0.2  V  increases  as  the  GDL  thickness  decreases,  and 


Fig.  6.  Effects  of  GDL  thickness  on  the  local  current  density. 
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decreases  along  the  axial  direction.  Nevertheless,  the  local 
current  density  for  V  =  0.7  V  remains  as  a  constant  and  do  not 
show  the  same  trend  as  that  for  V  =  0.2  V  as  mentioned  previ¬ 
ously.  It  should  be  noted  that  the  model  used  in  this  analysis 
is  isothermal  and  single  phase  condition,  the  temperature  gra¬ 
dient  effect  of  the  GDL  and  two-phase  condition  on  the  GDL 
shall  be  discussed  in  the  future  work.  Since  some  have  shown 
a  thinner  GDL  is  more  easily  flooded  due  to  lower  storage  of 
liquid,  thus  performance  under  2-phase  conditions  is  better 
for  a  thicker  GDL.  Also,  temperature  gradient  effects  of  the 
GDL  are  very  important  as  well  as  electron  transport. 

4.3.  Effect  of  flow  direction 

In  this  study,  the  effects  of  flow  direction  of  reactant 
fuels  in  both  the  anode  and  cathode  sides  on  the  cell  per¬ 
formance  are  also  investigated.  Fig.  7  represents  the  effects 
of  the  GDL  porosity  on  the  cell  performance  with  the  reac¬ 
tant  fuels  flowing  in  anode  and  cathode  sides  being  co-flow 
and  counter-flow  conditions.  It  is  obvious  from  Fig.  7  that 
the  cell  performance  increases  as  the  GDL  porosity  increases 
for  the  counter-flow  condition.  It  is  also  noted  that  the  cell 
performance  with  counter-flow  condition  is  better  at  higher 
operating  voltage  conditions,  while  the  cell  performance  with 
co-flow  condition  is  better  at  lower  operating  voltage  con¬ 
ditions.  This  is  due  to  the  remaining  fuel  reacts  with  air 
under  higher  operating  voltage  conditions  with  counter-flow 
condition,  which  results  in  a  better  cell  performance.  This 
phenomenon  can  be  confirmed  by  the  effects  of  flow  direc¬ 
tion  of  the  reactant  fuels  on  the  local  current  density.  Fig.  8 
presents  the  distributions  of  the  local  current  density  along 
the  axial  direction  with  counter-flow  condition  for  four  GDL 
porosities  at  V  =  0.7  and  0.2  V.  Comparison  of  the  correspond¬ 
ing  results  in  Fig.  4  and  Fig.  8  indicates  that  the  local  current 
density  with  counter-flow  condition  is  higher  than  that  with 
co-flow  condition  at  V=  0.7  V.  Besides,  the  local  current  den- 


Fig.  7.  Effects  of  flow  direction  of  reactant  fuels  in  both  the  anode  and 
cathode  sides  on  the  polarization  curves  of  the  PEMFC. 


Fig.  8.  Effects  of  GDL  porosity  on  the  local  current  density  with  counter¬ 
flow  condition  of  reactant  fuels  in  both  the  anode  and  cathode  sides. 

sity  with  counter-flow  condition  at  V=0.2V  is  higher  near 
the  entrance.  However,  the  local  current  density  with  co¬ 
flow  condition  becomes  higher  than  that  with  counter-flow 
condition  at  the  downstream  (F>  0.03  m).  Because  the  elec¬ 
trochemical  reaction  rate  is  high  under  lower  operating  volt¬ 
age  conditions,  the  fuel  consumption  rate  increases  near  the 
entrance  with  counter-flow  condition,  which  in  turn,  causes 
the  increase  in  local  current  density  and  less  fuel  along  the 
axial  direction. 


5.  Conclusion 

A  steady,  a  two-dimensional,  isothermal,  multi- species, 
full-cell  numerical  model  to  analyze  the  performance  of 
PEMFCs  has  been  established.  Liquid  water  effects  were  con¬ 
sidered  in  this  work.  The  effects  of  GDL  porosity,  membrane 
thickness,  GDL  thickness  and  flow  direction  of  reactant  fuels 
on  the  reactant  gas  transport  and  the  cell  performance  are 
examined  in  detail.  Brief  summaries  of  the  major  results  are 
listed  in  the  following: 

1.  The  cell  performance  is  enhanced  as  the  GDL  poros¬ 
ity  increases  with/without  consideration  of  liquid  water 
effect.  This  is  because  more  reactant  gases  transfer  into 
the  catalyst  layer  with  higher  GDL  porosities,  which  in 
turns,  more  chemical  reaction  occurs  and  more  reactant 
gases  consume,  and  a  better  performance  for  the  fuel  cell 
is  obtained. 

2.  The  cell  performance  increases  with  a  decrease  in  GDL 
thickness  at  lower  operating  voltage  conditions.  This 
is  due  to  the  fact  that  higher  concentration  gradients 
are  resulted  from  the  decrease  in  the  GDL  thickness, 
which  in  turns,  results  in  a  higher  mass  transfer  to  the 
catalyst  layer  and  then  more  electrochemical  reaction 
occurs. 

3.  The  cell  performance  with  counter-flow  condition  is 
higher  at  a  higher  operating  voltage  condition,  while  the 
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cell  performance  with  co-flow  condition  is  higher  at  a 
lower  operating  voltage  condition. 
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